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一 种 根据 羽 流 扩张 角 判 断 火 箭 发 动机 海平 面 排 气 速度 的 方法 


REM 


峰 “， 郭 常 超 


Et 电 集成 器 件 研究 所 ， 苏 州 215163; 
2. 中 国民 用 航空 飞行 学 院 飞行 技术 学 院 ， 广 汉 618307: 
3. 西安 近代 化 学 研究 所 ， 西 安 710065) 
摘要 : 近年 ,火箭 发 动机 的 非 接触 式 测量 研究 受到 各 国 的 关注 , 但 是 主要 关注 第 发 
动机 的 羽 流 湿度 场 或 振动 力学 研究 ,并 未 用 于 测 基 火箭 发 动机 的 排 气 速度 。 传统 排 气 速度 测 
量 方法 设备 复杂 , 且 长 期 没有 其 他 方法 作为 因此 本 文 提出 了 一 种 利用 羽 流 图 像 中 扩张 
角 测量 结果 ， 估 算 火 生发 动机 海平 面 排 气 速度 的 方法 ， 该 方法 结合 一 维 可 压缩 绝热 ; 


排 气 速度 和 
该 估算 方法 上 


火箭 发 动机 射流 : 海平 面 排 气 速度 ， 羽 流 阔 扩 ; 


比 冲 
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Abstract: Global researchers have recently become interested in non-contact measurements of 
rocket motors, but the primary focus is on the rocket motor's temperature field or mechanical 
vibration. As a result, non-contact methods were seldom used to measure the exhaust velocity of 
the rocket motor. Traditional exhaust velocity measurements need complex equipment and have 
had few references for a lengthy time. Thus, a method for estimating the rocket motor sea-level 
exhaust velocity was offered, which makes use of the plume expansion angle measured in its 


image by integrating a one-dimensional compressible adiabatic flow model, the Prandtl-Mayer 
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function in the ambient atmosphere, and the Brent algorithm. The relationship between the plume. 


expansion angle and the sea-level exhaust velocity was determined in this study, and the results 


were compared to those obtained using the traditional method to verify the new method's accuracy. 


The estimation results showed that there was a linear relation between the sea-level exhaust 


velocity and the plume expansion angle, and the relative error compared to the tradition method is 


less than 2%. Therefore, it was shown this estimation method has a fairly good accuracy. 


Key words: non-contact measurement; rocket motor jet flow; sea-level exhaust velocity; 


plume expansion angle; specific impulse. 
1 绪 言 
无 接触 式 测量 具有 不 影响 被 测 场 ,量程 上 限 
高 ,测量 延迟 低 等 优点 , 引起 许多 研究 者 的 关注 。 
但 是 在 火箭 发 动机 的 相关 研究 中 , 大 多 集中 在 火 
箭 发 动机 羽 流 温度 场 测量 或 振 WRA, 
航天 科 工 六 院 提出 一 种 凝 相 冉 粒 通过 该 
种 测量 方法 可 以 得 到 羽 流 当中 的 凝 相 物 颗粒 的 
运动 速度 59， 但 是 凝 相 颗 粒 速度 和 火箭 发 动机 
的 排 气 速度 并 不 相同 - 


而 传统 的 比 冲 或 排 气 速度 测量 方法 利用 推 
这 荐 计 和 称 重 台 等 大 型 复杂 设备 Cs1， 对 
进入 推力 室 前 的 火箭 燃料 进行 称 重 ， 


结束 后 得 到 平均 比 冲 或 者 平均 排 气 速度 -因此 


传统 方法 不 仅 操作 设备 和 过 程 复杂 、 原 理 单一 、 
无 参照 结果 , 并 且 发 动机 瞬时 排 气 


速度 。 0 

火箭 发 动机 羽 流 在 出 口 截面 存在 内 外 压 差 ， 
在 离开 喷 管 后 会 Tic, 因此 属于 一 种 超声 速 
流体 折 转 现象 - 综 上 , 结合 无 接触 式 测量 的 优点 ， 
为 丰富 火箭 发 动机 测量 手段 ,拓展 可 测量 物理 量 
范围 ,本 文 从 超声 速射 流 折 转角 和 马赫 数 关系 的 
角度 出 发 , 以 羽 流 图 像 为 依据 、 结合 出 口 截面 当 
地 声速 和 一 维 可 压缩 绝热 管 流 模型 ,提出 一 种 判 
断 海平 面 排 气 速度 的 估算 方法 .并 将 其 结果 与 传 
统 海平 面 排 气 速 度 测量 方法 的 结果 对 比 , 以 验证 


估算 结果 的 有 效 性 。 


2 判断 具体 方法 和 误差 来 源 
2.1 估算 理论 和 假设 


将 火箭 发 动机 的 羽 流 简化 为 连续 桶 型 波 串 
12, Anpa 1。 此 种 射流 形状 是 典型 的 低 度 欠 膨 胀 
射流 时 ( 即 出 口 截 耐 主流 压强 pe 和 环境 压强 pa 
之 比 1<n<1.15), 因 此 其 羽 流 主流 与 发 动机 中 轴 
线 在 第 一 个 波 串 单元 当中 有 一 个 初始 折 转 角 5- 
根据 普 朗 特 - 迈 耶 关 系 ， 此 转折 角 和 马赫 数 具 有 
一 一 对 应 关系 ， 所 以 通过 普 朗 特 - 迈 耶 函 数 的 逆 
解 能 够 得 出 出 口 截面 羽 流 的 马赫 数 M 121, 再 根 
据 当地 声速 we 就 可 以 得 到 羽 流 喷射 速度 上 入 
出 口 截面 温度 T. 的 函数 关系 。 同 时 ， 一 般 情况 
下 , 因为 火箭 发 动机 喷 管内 流动 速度 极 快 , 是 绝 
热 良好 的 冻结 流 ， 所 以 总 烩 守恒 上 0， 据 此 能 够 
建立 ve、7: 和 推力 室温 度 的 关系 。 综 上 ， 将 
EUR-A NGA FER, ARAA 
流 初始 折 转 角 5 求 出 海平 面 排 气 速度 we。 


Fig.1 Plume Wave Series 


itii ki E 
下 关系 pn 
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式 中 : PARRA, erg Jg: 
压 比 热 ,单位 为 J(kgK), 在 这 里 将 火箭 发 动机 燃 
气 视 作 量 热 理 想 气体 -所 以 对 于 同一 台 发 动机 而 
言 , 出口 截面 和 推力 室 之 间 的 能 量 关系 可 以 表示 
WF: 


[^] 


式 中 : 下 标 0 代表 推力 室 参 数 , 下 标 e 代表 出 口 
截面 参数 。 与 出 口 射 流 相 比 , 推力 室 的 燃气 
速度 很 低 ， 所 以 w 可 以 在 式 (2) 中 忽略 ， 则 式 
D 简化 及 变形 为 如 下 形式 : 


[3] 


对 于 理想 气体 而 言 定 压 比 热 m、 定 容 比 热 
cv、 绝热 系数 上 气体 常数 尺 的 关系 可 以 表示 如 
EN 


tak (4) 
R=R,/M 
式 中 : 通用 气体 常数 Re=-8.314J (mol * K), M 
为 该 气体 的 摩尔 质量 ， 单 位 为 kg/mol。 
结合 马赫 数 定义 M。 evi v, 和 当地 声速 的 
表达 式 w=VERF ， 式 (3》 和 式 (4) 联 立 可 以 
得 到 


G) 


uz) 


式 中 : Mos 即 为 出 口 截面 马赫 数 ，wao 为 推力 室 
的 当地 声速 。 根据 普 朗 特 - 迈 子 函数 ， 图 1 所 示 
的 右 伸 波 系 中 M6 和 出 口 截面 的 羽 流 扩散 角 5 
关系 可 以 具体 描述 如 下 : 


k+l 
sp Je Elem E a 


—arcian |M, 


一 般 认为 式 (6) 不 存在 普 适 的 连续 可 微 的 道 函 
KOEN, B M, mM. (OB TE CE E I A EU 
所 以 本 文采 用 布 伦 特 法 求 Mae Ms TE Be 


m 


2.2 误差 来 源 分 析 
(D 测量 误差 
加 取景 角度 和 测量 角度 


因为 摄像 者 和 火箭 发 射 台 角度 相对 位 置 关 
系 ,所 以 在 图 像 当中 无 法 准确 测量 羽 流 切 向 角度 
OUS. 


JL enn 
pon 
羽 流 边界 


仰视 拍摄 俯视 拍摄 
Fig2 Errors caused by inclined viewing 

回 羽 流 的 发 光 组 分 

对 于 液体 火箭 发 动机 而 言 ， 组 分 和 反应 


产物 透明 度 较 高 。 但 是 对 于 固体 火箭 发 动机 而 
流 当中 存在 着 铝 、 镁 等 物质 ， 过 程 中 
会 产生 峻 眼 的 亮光 , 所 以 固体 火箭 发 动机 的 羽 流 
eb d EL a 

[Sp ETR 

部 分 型 号 当中 存在 因为 发 动机 内 部 发 生 碍 
塞 导 致 的 羽 流 湾流 [3， 使 得 羽 流 结构 发 生 非 轴 
对 称 变化 ， 如 图 3 所 示 A-4 火箭 发 动机 
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同时 可 能 会 存在 羽 流 激 波 弓 [0， 改 变 喷 管 附近 
的 室 气 折射 度 ， 从 而 影响 测量 。 


Fig3 Asymmetric Vortexes near an A-4 Rocket 


Motor Nozzle Ol 
(2) 气体 参数 和 发 动机 设计 


本 文 所 涉及 的 气体 参数 燃烧 温度 Tn、 绝热 
系数 上 和 摩尔 质量 M 来 参考 自 [1 和 [23]. 的 
参考 值 是 在 一 定 压强 下 、 组 分 为 化 学 恰当 油气 比 
的 条 件 下 通过 最 小 自由 能 法 计算 得 来 Pr23, 由 于 
燃烧 室 工作 压强 和 初始 点 火 阶段 油气 比 与 理论 
假设 具有 一 定 偏差 ， 所 以 和 实际 火箭 发 动机 
的 推力 室温 度 并 不 相同 。 [I 火箭 发 动机 
喷 管 并 非 完全 绝热, ZE EUR ACE HEHE 
力 室 有 所 下 降 - 并 且 液 体 火箭 发 动机 推力 具有 可 
调节 性 ， 而 固体 火箭 发 动机 同类 推进 剂 组 分 多 
祥 ， 因 此 本 文 所 取 的 上 和 M 和 实际 工作 状态 有 
一 定 差距 。 ， 药 型 、 热 力 循 环 、 内 部 流 场 布 
局 等 设计 方案 较为 复杂 , 可 能 会 对 出 口 截面 流体 
性 质 产生 影响 。 


(3) 数值 求解 过 程 的 误差 


由 于 采用 了 布 伦 特 方法 对 普 朗 特 - 迈 耶 函 数 
REGERE, 所 以 产生 了 关于 Moe 的 后 向 误差, 为 
了 求 取 足 够 精确 的 出 口 截 面 马赫 数 , 本 文 的 后 向 
VOEUDGEIS EM.) 1.0X 105. 


2.3 求解 算法 

综 上 ， 海 平面 排 气 速度 估算 过 程 分 为 两 步 
外 利用 普 朗 特 - 近 耶 函数 逆 解 通过 布 伦 特 方法 求 
出 出 口 截面 马赫 数 ,加 利用 一 维 绝热 气体 动力 学 
关系 求 出 海平 面 排 气 速度 。 输 入 参数 分 为 三 类 : 
@ 气 动 几何 参数 ， 如 羽 流 扩张 角 o. 回 推力 室 气 
体 热力 学 参数 ， 如 推力 室 燃 气温 度 To、 绝 热 系 
数 上 和 摩尔 质量 Ms 加 数值 计算 参数 , 如 后 向 误 
X E. 以 上 所 述 可 见于 图 4- 


e» 


直面 高 度 火 入 发 | 
动机 羽 流 图 像 


pu 
peu FATRA 
5 
EURU 
[1:7 
CELLED] (mr 
Tak M 


EXC] 


Fig4 Estimation Progress of Sea-Level Exhaust 


Velocity 


3 估 测 结果 和 传统 方法 测量 结果 对 
比 


现 取 复 合 推进 剂 固体 火箭 发 动机 、 肝 类 液体 
火箭 发 动机 和 乙醇 液体 火箭 发 动机 各 四 个 型 号 ， 
其 羽 流 图 像 、 扩 张 角 和 型 号 分 别 可 见 图 5、 图 6 
和 图 7- 
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6=19° 


=16.6° 
a. UA1205 b，CZ-l1 第 一 级 
6=16.4° 
€. MGM3IB d. MSSAI 


Fig.5 Solid Composite Propellant Rocket Motors" 


Plumes and Expansion Angles 21 


|. 6710.5? D 
a YR2I b. DaFY6-2 
c RD-253 do LR-87-AJ-5 


Fig6 Liquid Hydrazine Propellant Rocket Motors's 


Plumes and Expansion Angles 27301 


A-7 


i 
n 


c. RD-100 d. RD-103 
Fig.7 Liquid Ethanol Propellant Rocket Motors's 


Plume and Expansion Angles P341 


根据 图 4, 估算 海平 面 排 气 速度 需要 推力 室 
燃气 热力 学 属性 参数 ， 见 于 表 1- 


Table 1 Gas Characteristics inside Thrust Chamber 1 


E 


燃烧 温 摩尔 质量 

bene BË TWK AARRE Mig mol* 
复合 推进 剂 。 337 147 2930 
HR 3300 124 2180 
高 浓度 乙醇 335 113 2288 


Tl S. 6. 7 所 示 的 火箭 发 动机 型 号 传统 测量 
方法 所 得 的 海平 面 排 气 速度 载 于 表 2。 
Table 2 Sca-Level Exhaust Velocities of Variable 


Rocket Motors Measured by Traditional Method. 


Bean 


型 号 海平 面 排 气 型 号 海平 面 排 气 
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men Bm! Table 3 Fractional Error between Plumc-Expansion. 
Method and Traditional Method 
UA1205 2315 — CZAUR—HD 2430 
型 号 相对 误差 型 号 相对 误差 
MGM3IB — 2145 MSSAL 2317 
UA1205 099% CZ-l1 第 一 级 。 1.86% 
YF-21 2540 DaFY6-2 2556 
MGM3IB 1.46% M55A1 0.63% 
RD-253 2798 — LR-87-AJ-5 2530 
YF-21 046%  DaFY62 186% 
A4 1959 A-7 2256 
RD-253 122% — LR-ST-ALS 0.08% 
RD-100 1990 RD-103 2160 
A4 0.19% A-7 0.39% 
RD-100 0.46% RD-103 12494 
根据 估算 算法 得 到 三 种 燃料 的 火箭 发 动机 
海平 面 排 气 速度 -扩张 角 关 系 曲线 和 表 2 所 述 的 
各 型 号 海平 面 排 气 速度 的 者 可 见 图 8. 图 8 显示 ”4 结论 
在 海平 结果 和 估 测 算法 结果 当中 ,海平 面 


排 气 速 度 和 羽 流 扩张 角 呈 线性 关系 .同时 根据 图 
8， 能 够 得 到 估 测 方法 和 传统 方法 结果 的 相对 误 
差 ， 见于 表 3。 从 表 3 可 以 得 出 ， 羽 流 扩张 角 法 
和 传统 测量 方法 所 得 结果 差别 不 大 , 相对 误差 在 
0.08%~1.86% 之 间 ， 因 此 相对 误差 小 于 2%， 据 
此 可 以 认为 羽 流 扩张 角 估 测 法 具有 较 高 的 精度 。 


10 15 
羽 流 扩张 角 /(”) 


Fig Plume vs, Sea-Level Exhaust Velocity 


Compared to Experimental Results 


本 文 将 推力 室 燃 气 特性 和 火箭 发 动机 射流 
特性 相 结 合 , 提出 一 种 根据 羽 流 扩张 角 判 断 火箭 
发 动机 海平 面 排 气 速度 的 估算 方法 .该 方法 包含 
了 一 维 可 压缩 绝热 流体 模型 、 大 气 背景 下 的 射流 
折 转 善 朗 特 - 迈 耶 函数 和 布 伦 特 算法 。 

以 复合 推进 剂 火箭 发 动机 、 腊 类 火箭 发 动机 
和 高 浓度 乙醇 火箭 发 动机 为 例 , 利用 该 估算 方法 
得 出 了 海平 面 排 气 速度 - 羽 流 扩张 角 曲 线 , 同时 ， 
通过 该 曲线 的 计算 结果 和 推理 架 测量 法 所 得 结 
果 对 比 。 经 过 以 上 研究 显示 : 


OD. 所 例 举 的 三 种 火箭 发 动机 的 海平 面 排 气 速 
度 和 羽 流 折 转 角 呈 线性 关系 ， 比 例 约 为 
66.6-87.03m/(s.^): 


(O0 估算 方法 和 传统 方法 相 比 ， 相 对 误差 不 高 

于 2%， 可 以 认为 该 方法 能 够 为 火箭 发 动机 海平 

面 排 气 速度 提供 一 种 精确 度 较 高 的 非 接触 测量 

的 理论 依据 。 
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